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Abstract 



In this paper we study an one-dimensional two-species exclusion model 
with open boundaries. The model consists of two types of particles moving 
in opposite directions on an open lattice. Two adjacent particles swap their 
positions with rate p and at the same time they can return to their initial 
positions with rate q if they belong to the different types. Using the Matrix 
Product Ansatz (MPA) formalism, we obtain the exact phase diagram of 
this model in restricted regions of its parameter space. It turns out that 
the model has two distinct phases in each region. We also obtain the exact 
expression for the current of particles in each phase. 
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I. INTRODUCTION 



The stationary state properties of one-dimensional driven diffusive systems are cur- 
rently of much research interest [1-5]. These systems exhibit very interesting cooperative 
phenomena such as boundary-induced phase transitions, spontaneous symmetry breaking 
and single-defect induced phase transitions which are absent in one-dimensional equilib- 
rium statistical mechanics. Many physical phenomena such as hopping conductivity, 
growth processes and traffic flows can also be explained by these models [6-9]. One of 
the most basic model is the Asymmetric Simple Exclusion Process (ASEP), which shows 
a rich behavior [10]. The ASEP is a model of particles diffusing on a lattice driven by 
an external field and with hard-core exclusion. Models with more than one kind of par- 
ticles have also been investigated. The ASEP in the presence of a second class particle 
(impurity) has provided a framework for the study of shocks (see [5] and [9]). Another 
model of this kind is the Partially Asymmetric Simple Exclusion Process (PASEP). In 
this model, particles are allowed to hop both to their immediate right and left sites with 
unequal rates. This model has been studied both with open boundaries [20] and in the 
presence of an impurity on a ring [11]. A multi-species ASEP has also been suggested, 
which seems to be a simple realization for real traffic [12]. 

In this paper we consider a model containing two types of particles on a lattice of length 
L with open boundary condition. The two types of particles, which we refer to them as " 
positive " and " negative " particles, move in opposite direction. The positive (negative) 
particles are injected (removed) from the left-most site of the lattice and are removed 
(injected) from the right-most site of the lattice. Every where through the lattice, two 
adjacent particels interchange their positions, unless they are both positive or negative. 
The system evolves according to an stochastic dynamical rule as follows. In each in- 
finitesimal time step dt the following events occur at each nearest-neighbour pair of sites 
1 (1 < i < L - 1): 



(+)(o) - 


- (0)(+) 


with 


rate 


1 


(o)H - 


- R(o) 


with 


rate 


1 


(+)(-)- 


> (-)(+) 


with 


rate 


P 


(-)(+) - 


- (+)(-) 


with 


rate 


1 



(1) 



where (+) and (— ) 

cates an empty site. Also, in each infinitesimal time step dt, the following events may 
occur at the first (i = 1) and the last (i = L) site of the lattice: 



At site 




with rate a 
with rate f3 



1 



At site i — L / (°) H ^ ra*e a f . 

A * Slt6 % ~ L \ (+) - (0) witfc rate /3 (2j 

For o = 0, this model reduces to the model introduced in [13,14] which using simula- 
tion data and doing exact calculations has been extinsively studied. These authors have 
shown that for certain values of the parameters a, (5 and p the symmetry of dynamics un- 
der interchange of positive and negative particles and of their directions is spontaneously 
broken. 

In reference [15] Alcaraz et al have studied the A-species stochastic models with open 
boundary and found their related algebras, which appear in the MPA formalism first 
introduced in [16]. Our model can be considered as a N = 2 case which will be studied 
in details. 

The process (1) has also been considered on a closed ring in [17-19]. It has been shown 
that when the density of positive particles is equal to the density of negative ones, de- 
pending on the values of the parameters of the model, three phases exist: a pure phase 
in which one has three pinned blocks of only positive, negative particels and vacancies 
(where the translational invariance is spontaneously broken); a mixed phase with a non- 
vanishing current of particles; and a disordered phase. Here we study the effects of the 
open boundaries. For certain cases {(5 — 1 or a — oo) we are able to solve our model 
exactly and find the modified phase diagrams (in comparison with q = case). We will 
show that for a = oo, where the system is devoid of vacancies, only two phases exist. 
In j3 — 1 limit the model has also two distinct phases in which the current of positive 
particles is equal to those of negative ones. 

This paper is organized as follows. In section 2 we will present the exact solution of the 
model for the case a = oo using the known results. We will also obtain the exact gen- 
erating function of the partition function of the model using the Matrix Product Ansatz 
(MPA) and calculate the current of the particles in /3 — 1 limit. In the last section we 
will compare our results with those obtained in [13] for q — 0. 



II. MATRIX PRODUCT SOLUTIONS 

In this section we will show that the stationary probability of the model defined in 
(1) and (2) can be obtained using the MPA for two specific cases a = oo and (3—1. 
According to the MPA formalism, the stationary probability P({C}) of any configuratin 
{C} can be written as a matrix element of a product of non-commuting operators. Be- 
fore reviewing this approach we define some notations. We introduce two occupation 
numbers, r« and 6i, for each site i, where = 1 if site i is occupied by a positive particle 
and Tj = otherwise. Similarly, 9i — 1 if site i is occupied by a negative particle and 
Qi = otherwise. Since the process is exclusive, so that each site of the lattice can only 
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be occupied at most by one particle, each configuration of the system is uniquely defined 
by the set of occupation numbers {r^, Now the normalized stationay state weight for 
a lattice of size L can be witten as: 



P({n, 9i}) = —(W\ HinD + M + (1 - Ti - 9i)E}\V). (3) 

The normalization factor Zl in the denominator of the equation (3), which plays a role 
analogous to the partition function in equilibrium statistical mechanics, is a fundumental 
quantity and can be calculated using the fact J2{ n ,e t } P{{ T n — 1- Thus one finds 

Zl = E ( W \ Iii^ D + M + C 1 - n ~ e i) E )\V) = {W\G L \V) (4) 

{TiA} »=i 

in which G = D + A + E. The operators D, A and E correspond to the presence of a 
positive, a negative particle, and a hole respectively. These operators with the vectors 
|V) and (W\ satisfy a certain algebra which will be discussed below. 



A. The limit a — > oo 

In this limit, as soon as a hole appears at a boundary site, it is removed. Therefore 
in the steady state the lattice will be empty of holes. Now the dynamical rules given by 
(1) and (2) reduce to 







(+)(-)->(-)(+) 


with 


rate 


V 






(-)(+)->(+)(-) 


with 


rate 


Q 


At 


site 


i = l (-)-(+) 


with 


rate 
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At 


site 


* = L (+)-.(-) 


with 


rate 


P 



Using the MPA one obtains the following quadratic algebra for this case 

pDA - qAD = D + A 

(3(W\A=(W\ (6) 
0D\V) = \V). 
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Now if one imagine the negative particles as holes, the problem reduces to the single- 
species PASEP with open boundaries and equal injection and extraction rates. As we 
mentioned, recently the PASEP has been studies widely with open boundaries [20]. Us- 
ing the results given there, we find two following phases in the thermodynamic limit 
(L -> oo): 

p-q 2 

The current of the positive particels J + is equal to the current of the negative ones 
J_ and has its maximum value 



U = J- = (7) 
Also the density of the positive particles < > in the bulk has a power law behavior 

The density of the negative particels < 0j > can be obtained using the equality 

< Ti > + < 9i >= 1. 

II. S-<\ 

p-q 2 

The current of the positive particels is again equal to the current of the negative ones 
and is given by 

J + = j_ = /3(i- JL). (9) 
p-q 

The density profile of the positive particles is linear in the bulk which is a consequence 
of the superposition of shocks 

<n>^^— + i{i-2^—). (io) 

p — q p — q 

This phenomenon has also been observed in the ASEP with open boundaries when the 
injection and extraction rates become equal [10]. 
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B. The limit (5 = 1 



Another limit which can be solved using the MPA formalism exactly is (3 — 1. The 
operators and vectors satisfy the following algebra 



Following [13] one can choose 



plJ j\ — qjiJJ 




DE 


= aE 


EA 


= aE 


E\V) 


= \V) 


D\V) 


= a\V) 


{W\E 


= (W\ 


(W\A 


= a(W\. 


-\V)(W\ , 


(V\W) = 1. 



Then the algebra (11) can be written as 



(11) 



(12) 



(-)DA — (-)AD = D + A 

a a 

D\V)=a\V) (13) 
(W\A = a(W\. 

This algebra is very similar to the algebra associated with the PASEP [19,20]. Here we 
adopt the same representation proposed in [19]. One can easily check that the following 
representation satisfy (12) and (13) 



D 



a 



p-q 



( 1 + a ^cl 

l + a(J) ^ 

l + a(f) 2 ^ 





l+«(f) 3 



, A = D T , E = \0)(0\. (14) 



V- • • ■ ■) 

Where the superscript T indicates the transpose, 
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c n = (l-( q -T)(l-a 2 ^) n - 1 ) , a = p-g-l, (15) 

(0| = [1000....] and |0) = (0| T . Using the matrix algebra given by (12) and (13) we find 
the following expressions for the current of the positive and the negative particles in the 
stationary state 

(W\G t '\pDA-qAD + DE)G L ' l - l \V) (W\G L ~ l \V) 
+ (W\G L \V) (W\G L \V) ' 

(W\G t ~ 1 (pDA-gAD + EA)G L - t - 1 \V} _ (WlG^Y) 
~ (W\G L \V) ~ a (W\G L \V) ' ( ] 

As can be seen from (16), the currents are site-independent (as it should be in the sta- 
tionary state), equal and given by the matrix element of powers of G. In what follows 
we will introduce a generating function to calculate the matrix element of all powers of 
G betwen the vectors \V) and {W\. 
Define a generating function 



/(^^A^WIV). (17) 

L=l 

The convergence radius of this formal series, R, is proportional to the current of particles 
given by (16) in the thermodynamic limit 

(W\G L - l \V) 

« = ""^~ V|ffl|v> • (18) 

On the other hand, the radius of convergence is the absolute value of the nearest singu- 
larity of /(A) to the origin. Once we obtain the singularities of the function /(A), we can 
calculate the current of particles and distinguish the phases. 
Using (12) and (13) one can expand the expression (VU|G L |U) as 



(G L ) := (W\G L \V) = ]T Yl (W\C ja E mi C jl ■ ■■E mr C jr \V) (19) 

r>0 . . . _ 

J0jr>0 

ji,mi,...,j, — i,m r >0 

jr'o+mi+jr'iH \-m r +jr=L 

where C := D + A. Noting that E m = E, after some computation, we obtain 
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(g l ) = (c L ) + E E 



[L-l-j -...-j r )\ 



r=l 

j0,jr>0 



1)!(L 



r ~ Jo 



3r)\ 



(C 30 ) . . . (C- 



'Jr\ 



31, 
jo+ji- 



■Jr-1>0 

—+jr<L-r 



(20) 



and also 



/(A) = E^ 1 (G i ) = 



L=l 



;i_ A -l) +A -2 (1 + A)E ^ =oA n + l (C n ) 

l-E~ A n+1 (C") 



(21) 



It is known that the expression g(X) := A n+1 (C n ) can be written explicitly in terms 
of the basic g-hypergeometric function (see [19] and references therein) 



g(x(X)) = x 



(p - q) \ P 



q r 2. q\ iq. q\ 



a 



(ax; 



"201 



ax, ax q q 

q ™2 ! _ 5 — 
-x p p 



(22) 



in which x(X) = \{ p -^ - 2 - V(^) 2 -4^}. The quantities (z;q) n and (z;^ are 
defined as 



(z;q)r 



fl, 2/ n = 0, 

\ (1-^)(1- ^(l-zg 2 ) ---(l- zq n ~ l ), if n>0, 

n=0 



(23) 



Lastly, the basic g-hypergeometric function is defined by the series 



2^1 



b 



E 

n=0 



(ai;g) w (q 2 ;g) r 
(6; g) n (g; g)„ 



(24) 



which tends to the usual hypergeometric series as q — > 1. The 20i series converges when 
< |g| < 1 and |z| < 1 [20]. In this paper the convergence condition of (22) is q < p; 
therefore, without losing the generality, we limit ourselves to this region. 
As we mentioned, the singularities of /(A) specify the phase diagram of the model. From 
the expression (21), we see that there are two possible sources for the singularities: the 
singularities of g(x(X)) and a zero of the denominator of (21). First we consider the sin- 
gularities of g(x(X)). From (22) one can see that g(x(X)) has two singularities: Ai = "^p- 



which is a square root singularity and A2 



p-g-j- 

a(p-q) 



which is a simple root. In order to 
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discuss the zeros of the denominator of /(A), we use the same assumption proposed in 
[19]. In the convergence region of (22) i.e. q < p, for < x < 1 the function g(x) 
satisfies g (x) > 0. It means that the function g(x) increases monotonically from to 
g{l) when < x < 1; therefore, the equation g{x{\^)) = 1 (which gives the zeros of the 
denominator of (21)) has only one root in this region. Comparing the absolute value of 
the singularities Ai, A 2 and A 3 , one can easily find the following results: 

/) For p — q > 2 the radius of convergence (18) of the formal series (17) is equal to 
R — A3 < Ai, A2 which is the solution of the equation g(x(\3)) = 1. Since A3 is a simple 
pole, we expect that Zl behaves asymptotically (L — > 00) as A3 L . The current of the 
particles (16) can also be obtained 



J + = J_ = a \ 3 (25) 

II) For p — q < 2 two different situations may occur. In the region specified by 
g(l) < 1 and p — q < 2, we find R — Ai = ^ and the current of particles to be 



J+ = J- = P -^ (26) 

In the region g{l) > 1, it turns out that R = X3 < X±, X2 which is again the solution 
of the equation (7(x(A 3 )) = I. The partition function Zl again behaves as and the 
current of particles in this case can be obtained from (25). The boundary of these two 
recent phases will be specified by 



9(1) = 



(p-q) 



a 



{p-q- I] 



q\2 



-2<t>i 



p — q — l,p — q — 1 q q 
- p p 



(27) 



In Fig.l we have plotted the phase diagram of our model in (3 — 1 limit both for q 7^ 
(the left diagram) and q = (the right diagram). As we mentioned in section (2) for 
a = 00, the line p — q > 2 is the line of shock configurations. The bold lines in Fig.l 
mark these lines. 
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FIGURES 
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FIG. 1. Plot of the phase diagrams for q = (the right diagram) and q ^ (the left 
diagram) in j3 = 1 limit. 

III. COMPARISON AND CONCLUDING REMARKS 

In this paper we studied a generalized two-species exclusion model with open bound- 
aries. The positive particles are supplied at the left end of the chain and they leave it 
at the right end. Similarly, the negative particles are supplied at the right end and they 
leave the system at the left end. As soon as a positive and a negative particle meet each 
other (the positive particle is supposed to be in the left hand side of the negative parti- 
cle), they interchange their positions with rate p. At the same time they may go to their 
initial positions with rate q. In q = limit this model reduces to the one studied in [13]. 
Using the MAP formalism, our model has been studied in two different limits (a = oo 
and (3 = 1) and the corresponding phase diagrams obtained. It has been shown that all 
the phases are symmetric in which the current of the positive and negative particles are 
equal. One can easily check that for q = all the results obtained here reduce to those 
obtained in [13]. In comparison with [13], as can be seen in Fig.l, the phase diagram of 
the model has been modified. 

The study of the whole parameters space of the model proposed in this paper is still 
an open problem. Using the mean field approximation and simulation data, the authors 
have shown that in q = limit this model has also asymmetric phases where the current 
of the positive and negative particles become different. It will be interesting to study the 
structure of asymmetric phases of this model which may occur for the certain values of 
the parameters a, (3, p and q. 
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